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The relativistic approach to electroweak properties of two- 
particle composite systems developed in Ref. [1] is generalized 
here to the case of nonzero spin. This approach is based on 
the use of the instant form of relativistic Hamiltonian dynam- 
ics. The generalization makes use of a special mathematical 
technique for the parametrization of matrix elements of elec- 
troweak current operators in terms of form factors. In this 
technique the parametrization is a realization of the Wigner- 
Eckart theorem on the Poincare group and form factors are 
reduced matrix elements. As in the case of zero spin the elec- 
troweak current matrix element satisfies the relativistic co- 
variance conditions and in the case of electromagnetic current 
it also automatically satisfies the conservation law. Physical 
approximations such as, for example, the relativistic impulse 
approximation, are formulated in terms of reduced matrix ele- 
ments. The electromagnetic structure of p meson is calculated 
as an example of realization of the technique proposed. 
PACS number(s): 13.40.-f, 11.30.Cp 

I. INTRODUCTION 

A new relativistic approach to electroweak properties 
of composite systems has been proposed in our recent 
paper [1] . The approach is based on the use of the instant 
form (IF) of relativistic Hamiltonian dynamics (RHD). 
The detailed description of RHD can be found in the 
review [2]. Some other references as well as some basic 
equations of RHD approach are given in Ref. [1]. 

In the paper [1] our approach was used to perform 
a realistic calculation of electroweak properties of pion 
considered as composite quark-antiquark system. The 
electromagnetic form factor and the lepton decay con- 
stant were calculated for pion using different model wave 
functions for the relative motion of quarks in pion. 

Now our aim is to generalize the approach to more 
complicated systems, namely, to composite systems of 
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two particles of spin 1/2 with nonzero values of total an- 
gular momentum, total orbital momentum and total spin. 
The main problem is a construction of electromagnetic 
current operator satisfying standard conditions (Lorentz 
covariance, conservation law etc., see, e.g., Ref. [1]). 

The basic point of our approach [1] to the construc- 
tion of the electromagnetic current operator is the gen- 
eral method of relativistic invariant parameterization of 
local operator matrix elements proposed as long ago as 
in 1963 by Cheshkov and Shirokov [3]. This canonical 
parametrization of local operators matrix elements was 
generalized to the case of composite systems of free par- 
ticles in Refs. [4,5]. 

In fact, this parametrization is a realization of the 
Wigner-Eckart theorem for the Poincare group and so it 
enables one for given matrix element of arbitrary tensor 
dimension to separate the reduced matrix elements (form 
factors) that are invariant under the Poincare group. The 
matrix element of a given operator is represented as a 
sum of terms, each one of them being a covariant part 
multiplied by an invariant part. In such a representation 
a covariant part describes transformation (geometrical) 
properties of the matrix element, while all the dynamical 
information on the transition is contained in the invariant 
part - reduced matrix elements. In the case of composite 
systems these form factors appearing through the canon- 
ical parameterization are to be considered in the sense of 
distributions, that is they are generalized instead of clas- 
sical functions. As was demonstrated in Ref. [1] this fact 
takes place even in nonrelativistic case. It is in terms of 
form factors that the electroweak properties of composite 
systems are described in the frame of the approach [1] . 

In our approach some rather general problems arising 
in the description of composite quark models have been 
solved. For example, the description of electromagnetic 
properties of composite systems in terms of form factors 
in Ref. [1], in fact, solves the problem of construction of 
the electromagnetic current satisfying the conditions of 
translation invariance, Lorentz covariance, conservation 
law, cluster separability and equality of the composite 
system charge to the sum of constituents charges (charge 
nonrenormalizability) . 
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Let us note that the importance of the problem of the 
construction of the electromagnetic current is actual not 
only for RHD but for all relativistic approaches to com- 
posite systems, including the field theoretical approaches 
[6-13]. 

We will construct the electromagnetic current opera- 
tor in the frame of IF of RHD. The same problem was 
considered in Rcfs. [6,13] in the point form of RHD and 
in Refs. [11,12] in the case of light-front dynamics. 

Our approach is a generalization of the method [13] to 
the case of the instant form dynamics. However, the sce- 
nario of the generalization of the Wigner-Eckart theorem 
is quite different. 

Physical approximations that we use in our approach 
are formulated in terms of reduced matrix elements, for 
example, the well known relativistic impulse approxima- 
tion. It means that the electromagnetic current of a com- 
posite system is a sum of one-particle currents of the 
constituents. It is worth emphasizing that in our method 
this approximation does not violate the standard condi- 
tions for the current listed before. To-day a construction 
of relativistic impulse approximation without breaking of 
relativistic covariance and current conservation law is a 
common trend of different approaches [6,10,11,13]. Let 
us note that in the present paper it is for the first time 
that such a construction has been realized for the case of 
nonzero spin in the frame of IF RHD. This is a variant 
of the relativistic impulse approximation (IA) formulated 
in terms of reduced matrix elements (see Ref. [1]) - the 
modified impulse approximation (MIA). 

The canonical parameterization of electroweak current 
matrix element, i.e. the extraction of the reduced matrix 
elements in the case of zero total spin and zero total an- 
gular momentum was performed in Ref. [1] and is rather 
simple. The case of composite systems with nonzero val- 
ues of total spin and angular momentum requires the 
development of a general method for canonical param- 
eterization of local operator matrix elements. Here we 
develop an adequate mathematics using as a base the 
paper [3]. 

In the present paper we propose a general formalism 
for the operators diagonal in the total angular momen- 
tum. The case of non-diagonal operators describing, for 
example, the radiative transitions between vector and 
scalar composite particles will be considered elsewhere. 
We demonstrate the application of the formalism in the 
case of a system with total spin one and total angular 
momentum one and with zero orbital momentum. In 
this connection the p - meson electromagnetic structure 
is calculated as an example. Using different model wave 
functions of the quark relative motion we calculate the 
electromagnetic form factors and the static properties of 
the p meson supposing quarks to be in the S state of rel- 



ative motion. It is interesting to mention that relativistic 
effects occur to produce a nonzero quadrupole momen- 
tum and quadrupole form factor. It is well known that 
in the nonrelativistic case the non-zero quadrupole form 
factor is caused by the presence of the D- wave and is 
zero otherwise. 

It is worth noticing that our approach guarantees the 
uniqueness of the solution for the p - meson electromag- 
netic form factors in contrast with the calculations based 
on the standard light-front dynamics in Refs. [14,15] 
where it is shown that the form factors obtained from 
matrix elements with different total angular momentum 
projections differ from one another themselves. An anal- 
ogous ambiguity takes place in the light-front dynam- 
ics calculations of the deuteron electromagnetic structure 
[16]. The next step of generalization of our approach con- 
cerns two-particle systems with nonzero orbital momen- 
tum giving a possibility to describe this simplest nuclear 
composite system [17]. 

The paper is organized as follows. In Sec. II the canon- 
ical parametrization of local operator matrix elements 
between one-particle states of arbitrary nonzero spin is 
described. This parametrization presents the extraction 
of the reduced matrix clement on the Poincare group. 
The electromagnetic current matrix element is derived 
for the spin 1/2 particle. For this case the relations be- 
tween the form factors in the canonical parametrization 
form and commonly used Sachs form factors are given 
explicitly. 

In Sec. Ill we show how to construct the electromag- 
netic current operator for composite system of two free 
particles of spin 1/2. The current matrix element is ob- 
tained in the basis where the center-of-mass motion is 
separated. The electromagnetic properties of the system 
are defined by reduced matrix elements, or the so called 
free two-particle form factors, these form factors being 
generalized functions. This means that, for example, the 
static properties of the system are given by the weak 
limits as Q 2 — > (Q 2 = —q 2 , q is the momentum trans- 
fer). A special attention is paid to the case of total spin 
one and zero total orbital momentum. In this case the 
electromagnetic properties are defined by four free two- 
particle form factors - charge, quadrupole, magnetic and 
the magnetic quadrupole form factor of the second kind 
(see the review in Ref. [18] and the references therein). 

In Sec. IV the electromagnetic current matrix element 
is constructed for the case of two interacting particles. 
Each step of the calculation in the developed method re- 
mains Lorentz covariant and conserves the current. The 
composite system form factors are derived as some inte- 
gral representations in terms of the wave functions ob- 
tained in the frame of RHD. 

In Sec.V the developed formalism is used in the case 
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of the system with total spin one and total angular mo- 
mentum one. In terms of reduced matrix elements the 
so called modified impulse approximation (MIA) [1] is 
formulated. In MIA p - meson form factors are obtained 
explicitly. 

In Sec. VI the results of calculations of static proper- 
ties and electromagnetic form factors of p meson are dis- 
cussed. 

In Sec. VII the conclusion is given. 

II. PARAMETRIZATION OF ONE-PARTICLE 
OPERATOR MATRIX ELEMENTS 

A parametrization of the current operator matrix ele- 
ment is a representation of the matrix element in terms 
of form factors. Let us remark that it is just form factors 
- Lorentz invariant functions - are extracted from scat- 
tering data. In simple cases, such as particles with spin 
or 1/2 the parametrization can be obtained through 
semiempirical approaches. However, in more complicated 
cases, for example, in the case of nuclei with arbitrary 
spins one needs a general mathematical method. 

Let us describe now a general method for canonical 
parameterization of local operator matrix elements ( sec 
Rcf. [3], too). As we have mentioned before, here we 
are dealing with the matrix elements diagonal in total 
angular momentum only. Non-diagonal case will be con- 
sidered elsewhere. 

From the group theory point of view a parameteriza- 
tion is, in fact, a realization of the well known Wigncr 
Eckart theorem on the Poincare group. The parame- 
terization extracts the reduced matrix elements that are 
invariant under the Poincare group. These reduced ma- 
trix elements are just form factors that are measured in 
experiments. 

The main idea of the parametrization can be formu- 
lated as follows. Using the variables entering the state 
vectors which define the matrix elements one has to con- 
struct two types of objects. 

1. A set of linearly independent matrices which are 
Lorentz scalars (scalars or pseudoscalars). This set de- 
scribes transition matrix elements non-diagonal in spin 
projections in the initial and finite states, as well as the 
properties defined by the discrete space-time transfor- 
mations. 

2. A set of linearly independent objects with the same 
tensor dimension as the operator under consideration (for 
example, four-vector, or four-tensor of some rank). This 
set describes the matrix element behaviors under the ac- 
tion of Lorentz group transformations. In the case of 
Lorentz scalar operator the second set coincides with the 
first one. 



The operator matrix element is written as a sum of 
all possible objects of the first type multiplied by all the 
possible objects of the second type. The coefficients in 
this representation as a sum are just the reduced matrix 
elements - form factors. The obtained representation is 
then modified with the use of additional conditions for 
the operator, such as the conservation laws, for example. 
In order to satisfy these additional conditions in some 
cases some of the coefficients occur to be zero. 

To demonstrate this let us consider the parameteriza- 
tion of the matrix elements taken between the states of 
a free particle of mass M in different simple cases. Let 
us normalize the state vectors as follows [1]: 

(p,m\p' ,m') = 2po5(p- p')5 mm > , (1) 

Here p , p' are three-momenta of particle, po = 
\J M 2 + p 2 , m , m! are spin projections. 

Let us consider first the parameterization of the matrix 
element of scalar operator A{x) taken between the states 
of free particle with zero spin. Because of translation 
invariance it is sufficient to consider A(0). Using the 
variables of state vectors in the initial and final states 
one can construct two linearly independent scalars: p 2 = 
p' 2 = M 2 , (p-p') 2 = q 2 = -Q 2 . Only the second one - 
the momentum transfer square - is nontrivial. Ignoring 
the trivial dependence on M we can write: 

(p,M\A(0)\p',M)=f(Q 2 ). (2) 

Let us consider now the matrix element of a scalar 
operator between the states of a particle with spin j: 

(p,M,j,m\A(0)\p',M,j,m') (3) 

The tensor dimension of the operator is not changed, 
nevertheless, the matrix element is now a matrix in spin 
projections in the initial and finite states. 

Let us construct the set of Lorentz scalars - linearly 
independent matrices in spin projections to be used for 
the construction of the representation of Eq. (3). Let 
us use the covariant spin operator T fJ, (p) [19]. In the 
rest frame this operator coincides with the particle spin 
operator: 

r°(0)=0, f(0) = Mj, [ji,jk]=ie ik iji. (4) 

The covariant spin operator can be defined with the use 
of the Pauli-Lubanski vector [20] . In terms of matrix 
elements we have: 

r w(p)Hp>™KlpX} 

= (0,m|J7- 1 (A p )w; A1 C/(Ap)|0,m') (5) 
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= (A p )^(0,mK|0,m') = (A P )r 1^,(0) . 

Here A p is the boost for the transformation from the rest 
frame of the particle to the laboratory frame, U(A p ) is 
the corresponding representation operator. The matrix 
A p can be written explicitly using the matrix for Lorentz 
transformation of the vector p' into the vector p: 



is 



M 2 



M 2 + p x p' x 



In our case A p is given by Eq. (6) with p' = 
(M, 0,0,0). 

Now equations (4) can be written as: 



r° mm ,(o) = (o,m\w \o,m') = o, 

f m m'(0) = {0 ,m\w\0 ,rri) = M j mm > . 
Using the explicit form of A p we obtain: 

P(pj) 



(7) 



r (p) = (p&'), r(p) = Mj + 



Po + M 

T 2 = -M 2 j(j + 1). (8) 

As is known [20], spin transforms under the action of 
Lorentz group following the so called little group which 
is isomorphic to the rotation group, that is the corre- 
sponding transformations are realized by the matrices of 
three dimensional rotations. The derivation of the ex- 
plicit form of these matrices with an arbitrary spin is 
described in Ref. [21]. For spin 1 and 1/2 the matrices 
have the form: 

D 1/2 ( Pl ,p 2 ) = cos^-2i(k-j) sin^ , 
D 1 (pi,p2) = I — i(k ■ j) sinw + (k ■ j) 2 (coslu — 1) , 



uj = 2 arctan 



[Pi ,P2_ 



(pio + Afi)(p2o + M 2 ) - {piPi) 

[Pi ,P2] 



k 



[Pi ,P2\ 



(9) 



So, the 4-spin operator is transformed under Lorentz 
transformations p^ = A^p' v in the following way: 



r<» = A<t^(p, p')T»(p')Dj(p',p) 



(10) 



Using Eq. (10) one can show directly that the ma- 
trix elements of the operators D^(p, p')T^{p') and 



F Ai (p)D J (p, p') transform as 4-pseudovectors, the ma- 
trix elements of the operators D^{p, p')p^^{p') and 
p' IJ T^(p)D : > (p, p') - as 4-pseudoscalars. 

We construct the set of linearly independent Lorentz- 
scalar matrices using the vectors p" 1 and the pseu- 
dovector D^{p, p')T^(p'). Note, that the pseudovector 
T^(p)D J (p, p') does not enter the decomposition, being 
linearly dependent. One can show this fact using the 
relation (10) and the explicit form Eq. (6) of A^, trans- 
forming p' into p . It is easy to obtain: 

T^p)D^p, p')=Di(p, p') [PV) 

p>* + p ip - 



M 2 + p^pl' 



\p^(p')} 



(11) 



As p' r M (p') = 0, the set in question of linearly indepen- 
dent matrices in spin projections of the initial and the 
final states giving the set of independent Lorentz scalars 
is presented by 2j + 1 quantities 



Di(p,p')(p^(p')) n , n = 0,l,...,2j 



(12) 



The number of linearly independent scalars in Eq. (12) 
is limited by the fact that the product of more than 2j 
elements r M (p') reduces as it is linearly dependent. If 
n is even, then the obtained quantities are scalars, if n 
is odd then they are pseudoscalars. The current matrix 
element (3) is represented by the linear combination of 
these linearly independent Lorenrz scalars. The coeffi- 
cients in this combination f n (Q 2 ) are just form factors. 
These form factors are invariant under rotations and so 
they do not depend on spin projections. So, they de- 
pend upon only one scalar combination of variables - the 
momentum-transfer squared. 

For self-adjoint operator A(0) a minor modification of 
the set (12) is necessary: in the scalar product in Eq. (12) 
the factor i appears. So, now the current matrix element 
(3) can be written in the form: 

(p,M,j,m\A(0)\p',M,j, to') 



2j j 

= E E (rn\D\p,p')\m") 

n—0 m" — —j 

x(m"\{ip^(p')r\m')f n {Q 2 ) (13) 

The obtained operator is self-adjoint. This can be easily 
shown using the following relation 

p'^(p)Di(p, p') = -Di(p, p')p^(p') (14) 

which is a consequence of Eq. (11). 
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In the case of a scalar operator the values n will be 
even and for a pseudoscalar they will be odd. 

Let us consider now the 4- vector operator ju(0). To 
parametrize the matrix element one needs a set of quan- 
tities of the appropriate tensor dimension. Using the vari- 
ables entering the particle state vectors one can construct 
one pseudovector T^(p') and three independent vectors: 



K /1 = {p- P% = % > K' fl = {p + p% , 



(15) 



Here e M „ A p is a completely anti-symmetric pseudo-tensor 
in four dimensional space-time with e 012 3 = —1. 

The operator matrix element contains the matrix el- 
ements of the listed quantities multiplied by (p, p') 
from the left. Each of such products is to be multiplied 
by the sum of linearly independent scalars constructed 
while obtaining the parameterization (13): 

(p, M,j,m\j^(0)\p', M,j, m') 



Y,(rn\D j (p, p')\m")(m"\ F 1 K>^ + F 2 r"(p') 

m" 

+F 3 R^ + F 4 K^\m') , (16) 



where 



2j 



(17) 



n=0 



Let us impose some additional conditions on the oper- 
ator. 

1. Let us require the operator to be self-adjoint. One 
can check this condition by use of Eq. (14). In the right- 
hand side (r.h.s.) of Eq. (16) we need to modify slightly 
the vector multiplied by F 2 . The new vector is a linear 
combination of the 4-vectors entering Eq. (16) and has 
the following form: 



FV) r<V) - ^72- WV)) 



(18) 



The terms containing F 2 and F 3 are modified in the fol- 
lowing way: 

FiA" -> ^(F i A» + A»F i ) = {F i A»} + , i = 2,3. 

(19) 

Here A^ are the vectors entering Eq. (16) and changed 
following Eq. (18). The terms containing F 3 and F4 have 
to be multiplied by i. 



2. It is useful to modify the parametrization so as to 
make the vectors - multipliers of Ft - orthogonal to each 
other. This results in the following form of the vector 
multiplied by F 2 (18): 



/ K' v K v \ 



)) 



(20) 



3. Let us impose the condition of parity conservation. 
This condition is satisfied if all the terms in the sum (17) 
contain an even number of pseudovector factors T^. So, 
the sums in F\ and F4 are over even n : 2j > n > 0, 
in F 3 over even n : 2j — 1 > n > 0, and in F 2 over 
odd n : 2j - 1 > n > 0. 

4. Let us impose the conservation condition j^K^ = 
j^q 11 = 0. It is easy to see that this condition is satisfied 
only if F 4 = 0. 

So, finally we have 

(p, M, j, m \jn(0)\p', M,j, to') 
= £ (m\D*{p, p')\m"){m"\ F X K^ + {F 2 [T»(p>) 

m" 

-W(P')) + f|)] } + + i {F 3 R,} + \m'), (21) 

This construction can be used, for example, to obtain 
the electromagnetic current matrix element in the case 
of particle with spin 1/2. 

Now let us list the conditions for the electromagnetic 
current operator to be fulfilled in relativistic case (see 
Ref. [1] and the references therein). 

(i).Lorentz covariance: 



(22) 



Here A is a Lorentz-transformation matrix, U(A) is an 
operator of the unitary representation of the Lorentz 
group. 

(ii). Invariance under translation: 



U- 1 {a)f{x)U{a)^f{x-a). 



(23) 



Here U (a) is an operator of the unitary representation of 

the translation group. 

(in) . Current conservation law: 



|4/(0)]=0. 



(24) 



In terms of matrix elements ( j^(0) ) the conservation law 
can be written in the following form: 



(25) 
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Here g M is four-vector of the momentum transfer. 

(iv) . Current-operator transformations under space-time 

reflections: 



= (j°{x°,-g),-'j{x°,-ar)) 
U R f{x)U^=f{-x). 



(26) 



In Eq. (26) Up is the unitary operator for the repre- 
sentation of space reflections and Up is the anti-unitary 
operator of the representation of space-time reflections 
R = PT. 

Our parametrization satisfies all the conditions listed 
above. The parametrization (21) for the electromagnetic 
current in the case of particle with spin 1/2 has the fol- 
lowing form: 

(p, M, 1 m|i M (0)|p', M, i, m!) 



= J2{m\D 1 / 2 (p, p')\m") 

m" 

x(m"\f 10 (Q 2 )K^+ t f 3Q (Q 2 )R^\m') , (27) 

The form factors f w (Q 2 ) = fi(Q 2 ) and / 3 o(Q 2 ) = 
f2(Q 2 ) are the electric and the magnetic form factors 
of the particle, respectively. These form factors can be 
rewritten in terms of standard Sachs form factors Ge{Q 2 ) 
(electric) and Gm{Q 2 ) (magnetic) [22]: 



h(Q 2 ) = 
/ 2 (0 2 ) = - 



2M 
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MyjAM 2 + Q 2 



G E {Q 2 ) 
G m (Q 2 ) 



(28) 



In analogous way (although a little more cumbersome) 
one can obtain the matrix elements of the operators of 
higher tensor dimension. This interesting problem, how- 
ever, is out of scope of the present paper. 



III. PARAMETRIZATION OF MATRIX 
ELEMENTS OF THE TWO-PARTICLE 
ELECTROMAGNETIC CURRENT OPERATOR 

To describe the properties of the system of interacting 
constituents in our approach it is necessary to have the 



reduced matrix elements on Poincare group (the form 
factors) which describe the properties of the composite 
system of free constituents. In this section we generalize 
the method of parametrization of the previous section to 
the case of such free systems (see Refs. [4,5], too). 

Let us consider a system of two free particles with spins 
1/2 and let us parametrize the matrix element describ- 
ing the transitions in this system. Let us construct, for 
example, the electromagnetic current operator matrix el- 
ement. The matrix element can be taken between the 
following two-particle state vectors: 



I Pi , mi; p 2 , m 2 ) =\pi, mi) <g> | pi , m 2 ) 



(29) 



Here pi , p 2 are three-momenta of particles, mi , m 2 
are spin projections on the z axis. The one-particle state 
vectors are normalized by Eq. (1). 

As well as the basis (29) one can choose another set of 
two-particle state vectors where the motion of the two- 
particle center of mass is separated: 

| P, y/s, J, I, S, mj) , 

(P, V7s, J, I, S, mj\P', V7, J', I', S', mj.) 

= N CG SW(P - P')5(y/i- Vs 1 )6jj>6wSss'6 mjmjl , 

(30) 



CG 



(2P 



8k* 



k 



y/A(s , M 2 , M 2 ) 
2v^ 



Here P M = (pi +^2)^, P 2 = s, y/s is the invariant mass of 
the two-particle system, I is the orbital angular momen- 
tum in the center-of-mass frame (cm.), S is the total 
spin in the cm., J is the total angular momentum with 
the projection mj; Mi , M 2 are the constituent masses, 
and A(a, b, c) = a 2 + b 2 + c 2 - 2(ab + be + ac). 

The basis (30) is connected with the basis (29) through 
the Clebsh-Gordan decomposition for the Poincare 
group. Here we write the decomposition in a little more 
general form than in Ref. [1]: 

\P, y/s, J, I, S, raj) 



El dpi 
, J 2pn 



Pi ,mi; P2 ,m 2 ) 



10 2p 2 o 



x(pi ,mi; p 2 ,m 2 | P, y/s, J, I, S, mj) , (31) 



where 



(pi,m 1 ;p 2 ,m 2 \P, yfs, J, I, S, mj) 
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= V^[X(s, Ml Mi)]- 1 ' 2 2P S(P - Pl - p 2 ) 
x 2<mi| D 1/2 (pi , p) | mi )(m 2 | D 1 ' 2 (p 2 , P) |m 2 ) 



x(- - Aims l-Sms) Yi m ,(i9,v?) (Slm s mi \ Jmj) . 

Here $ , <p are the spherical angles of the vector p = 
(pi — p 2 )/2 in the cm., YJ m; <p) are spherical har- 
monics, (1/2 1/2 rhi m2 \Sms ) , (S lm s mi\Jmj) are 
Clebsh-Gordan coefficients for the group SU(2), arc 
the known rotation matrices to be used for correct rela- 
tivists invariant spin addition [5] , the summation is over 
mi,m,2,mi,ms- 

The decomposition in spherical harmonics and angular 
momenta summation in Eq. (31) are performed in the 
cm. and the result is shifted to an arbitrary frame by 
use of D-functions [5] . 

The electromagnetic current matrix element for the 
system of two free particles taken in the basis (29) can 
be written as a sum of the one-particle current matrix 
elements: 

(pi , mi ; p 2 , m 2 (0) \p[ , mi ; p' 2 , m' 2 ) 

= (P2,m 2 |p 2 ,m 2 )(pi,mi|ii, J (0)|pi,mi) + (1 <-> 2) . 

(32) 

Each of the one-particle current matrix elements in 
Eq. (32) can be written in terms of form factors as in 
Section II. In the case of the particles with spin 1/2 wc 
make use of Eq.(27). So, in this case the electromagnetic 
properties of the system are defined by the form factors 
fi , h, given by Eqs. (27), (28). 

Now let us construct the electromagnetic current ma- 
trix element for the system of two free particles in the 
basis (30) following the previous Section. Let us con- 
sider first a simple case J = J' = S = S' = I = I 1 = 0. 
(We omit these variables in the state vectors.) This set 
of quantum numbers appears, for example, in the case of 
pion. Now there is no pseudovector r M , but along with 
the scalar (P — P') 2 = — Q 2 two additional nontrivial 
scalars do appear s' = P' 2 and s — P 2 - the invari- 
ant mass squares for the free two-particle system in the 
initial and in the final states. So, the form factors en- 
tering the parametrization are functions of the variables 
Q 2 , s , s' . The current matrix element is presented by 
the linear combination of the four-vectors P M and P^: 

(P,VS\jj?\o)\P',^) 

= {P» + Pl) 9i(s, Q 2 ,s') + {P„ - Pl) g 2 (s, Q 2 ,s') . (33) 



Making use of the conservation condition (25) 

j<°>(0)(P-Py = 0. (34) 
we can write the parametrization in the form: 

{P,Vi\jj? ) (0)\P',V?) 

= A^s,Q 2 ,s')g (s,Q 2 ,s'). (35) 

Here <?o(s, Q 2 , s') is the reduced matrix element. We will 
refer to this invariant as to the free two-particle form 
factor. The vector A ll (s,Q 2 , s') is defined by the cur- 
rent transformation properties (the Lorentz covariance 
and the conservation law): 

A^s,Q 2 ,s') 

= (1/Q 2 )[( S -s' + Q 2 )P M + (s' - s + Q 2 )P'J . (36) 

The free two-particle form factor can be expressed in 
terms of the one-particle form factors (27) (see Rcf. [1] 
for details). To do this one has to perform in Eq. (35) 
the Clebsh-Gordan decomposition of the irreducible rep- 
resentation (30) into the direct product of two irre- 
ducible representations (29), (31) and to take into ac- 
count Eqs. (27), (32). As the form factors are invariants 
one can perform the integration in Eq. (35) in the coor- 
dinate frame with P' = , P = (0, 0, P). The explicit 
form of go{s, Q 2 , s') in the case of two particles with spin 
1/2 and mass M can be found in Ref. [1]. 

Let us perform the analogous parameterization for the 
set of quantum numbers in the basis (30) in the case 
J , J' ^ (see Ref. [5], too). Finally we will take J = J'. 

The Lorentz covariant properties of the matrix clement 
are defined (in analogy with Eq. (15)) by three 4-vectors 
and one pseudovector: 

Kl = (P + P')^ i^ = (P-P%, 

Rn = e^ Xp p»p' x rr(p') , ryp'). (37) 

The pseudovector T^(P) does not enter the parametriza- 
tion because it can be expressed through r M (P') by the 
equation analogous to Eq. (11): 

-1= r(P)„ d*{p, p') = &(P, p') f -^W) 

V s I Vs' 

-^.^±45[ W po]1. 

^ p v p'» + VU j 

The set of linearly independent matrices to be used for 
the decomposition of the current matrix element is ob- 
tained from the vectors P M and T^(P') following Eq. (12): 
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D J (P, P')(P^(P')) n , n = 0,1,..., 2.7. (38) 

Using the conditions of self-adjointness, current con- 
servation (34), parity conservation and orthogonality 
of the parametrization vectors to one another (as in 
Eq. (21)) we obtain (see Eq. (19), too): 

(P, y/i, J, I, S, mj\jW\P', V?, J, I', S', m'j) 



= J2(m.,\D J (P,P')\m['j} 

m'J 

:(™:;iEK ss '^(^ 2 . s ')} lm ' jh (39) 
»=i + 



4 = Q2 [(« - s ' + Q 2 ) p » + - s + Q 2 ) p 'A . 



-(V5 + V7)^ 



+ 



Vg 7 -^ + yip; 
pp' + V^ 7 



A(s,-QV) 



(V^+ V^) 2 + Q 5 



Pp • 



[PTV)] , 



(40) 



The quantities Ff ss in Eq. (39) are defined by the re- 
lations analogous to Eq. (17): 



2J 

F? SS ' =J2 ft SS '{s, Q\ s')(iP^(P')r . (41) 

n=0 

The sum in Eq. (41) is taken using the parity conservation 
condition as in Eq. (21). 

Let us remark that the reduced matrix elements - the 
invariant form factors — now (in contrast with the form 
factor (35)) depend on the additional invariant quantities 
I, I', S, S' that are invariant degeneration parameters in 
the basis (30). 

The self-adjointness condition is fulfilled in Eq. (39) if 



f ll ' SS l(s,Q 2 ,s') = fir S (s',Q 2 ,s). 
Here the star means the complex conjugation. 



(42) 



Let us consider especially the case Mi = M 2 = 
M , J = J' = S = S' =1 , I = V =0. We will use the 
reduced matrix elements for this case to calculate p 
meson properties in MIA neglecting the D - state con- 
tribution (see, e.g. Ref. [15]). In this case the functions 
pii SS j n j-^-j have the following form (compare to 
Eqs. (17), (21)): 

Pi = f W (s,Q 2 ,s') + f 12 (s,Q 2 ,s , )(iP v V(P')) 2 , 

F 2 = .f 2 i(s,QV) (iP^(P')) , P 3 - / 3 o(s,QV) • 

(43) 

In equations (43) the fixed variables l,l',S,S' are omit- 
ted. Time reflection invariance imposes the following 
conditions: 

fL = fin, i = l,3; / 2 *i = -/ 2 i. (44) 

With Eqs. (42) and (44) we obtain 

fin(s, Q 2 , s') = fin(s', Q 2 ,s) , i = 1, 3 ; 

f 2 i(s,Q 2 ,s') = -f 21 (s',Q 2 ,s) . (45) 

The relation (45) demonstrates that the form factor /21 
appears in the parametrization as a consequence of the 
fact that the invariant masses in the initial and the final 
states are different. This form factor gives no contri- 
bution to clastic processes, for example, to the electron 
scattering by the composite particle, however, it docs 
contribute to radiative transitions. 

Let us rewrite Eq. (43) in terms of standard form fac- 
tors instead of /, n : 

Pi = .goc(s,QV) + SoQ(s,QV){(*P^(P')) 2 - 



-(l/3)Sp(W(P')) 2 } 



S P (P,L-(P')) 2 



P 2 = <7omq(s,QV) (iP.r^(P')) , 

P 3 = 9om{s,Q 2 ,s') . (46) 

The scalar factor in the term with <?oq is chosen in 
such a way as to make it possible to interpret goQ as 
a quadrupole form factor of the system of two free par- 
ticles. In other terms goc is the charge form factor, 
9omq — /21 is the magnetic quadrupole form factor of 
the second kind, its classical analog being the so called 
toroidal magnetic moment [18], goM = /30 is the mag- 
netic form factor of the free two- particle system. 
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In the same way as in the case of Eq. (35) (see Ref. [1] 
for the details) the free two-particle form factors in 
Eq. (46) can be written in terms of the constituent form 
factors. The corresponding equations are rather compli- 
cated, so we present here only the equations that we will 
use for the p-meson electromagnetic structure, that is for 
the case when the constituents are the u- and d- quarks 
(see Appendix 1). 

The free two-particle form factors in Eqs. (35), (39), 
(43) , and (46) are to be considered in the sense of distri- 
butions [1]. For example, g (s , Q 2 , s') in Eq. (35) has to 
be interpreted as a Lorentz invariant regular generalized 
function on the space of test functions 5(R 2 ) [23]. 

Let us define the functional giving the regular general- 
ized function as 

(go(s,Q 2 ,s') ,tp(s,s')) 

= J dn{s,s')g {s,Q 2 ,s')tp{s,s') . (47) 

Here 

dfi{s, s') = 16 9(s - AM 2 ) 9(s' - AM 2 ) 

x\fss' d[i(s) dn(s') , d^i(s) = ^kd^/s . (48) 

The quantity Q 2 is a parameter of the generalized func- 
tion, Mi = M 2 = M. The function 6(x) is the step 
function. 

ip(s , s') is a function from the space of test functions. 
So, for example, the limit as Q 2 — > (the static limit) 
in g (s , Q 2 , s') (35) exists only in the weak sense as the 
limit of the functional: 

lim (90,^) = ((e q + eg)S{ji(s') - fi(s)),ip) . (49) 

e q and e 5 are the constituent charges, 5 is the Dirac delta- 
function. 

It is just the limit in the sense (49) that gives the elec- 
tric charge of the free two-particle system. The ordi- 
nary point -wise limit of the form factor go &s Q 2 — >0 is 
zero. The equations analogous to Eq. (49) are valid for 
the static limits of the free two-particle form factors in 
Eqs. (46), (Al), (A2) (A3), too. 

Let us note that the conditions imposed on the free 
two-particle form factors that follow from the conditions 
for the electromagnetic current operator (42), and (45) 
have to be considered in the weak sense, too. 



IV. PARAMETERIZATION OF THE CURRENT 

OPERATOR MATRIX ELEMENTS FOR 
SYSTEMS OF TWO INTERACTING PARTICLES 

In this Section we generalize the parameterization 
method of the previous sections to the case of composite 
system with the structure defined by the interaction of 
its constituents. 

Let us consider the operator j p (0) that describes a 
transition between two states of a composite two- con- 
stituent system, ,7^(0) being diagonal in the total angular 
momentum. Let us neglect temporarily the additional 
conditions of self-adjointness, parity conservation etc. in 
the same way as when constructing the matrix elements 
(16) in Section II. The Wigner-Eckart decomposition of 
the matrix clement has the form (16), (17). To empha- 
size the fact that the particle is composite, let us rewrite 
Eqs. (16), (17) using new notations: 

(Pc m Jc \ j ^(O)] p c ',m' Jc ) 

= ( m Jc \D J c (p C) p' c ) [ F? Kl + F 2 C r»(p' c ) 

+FiR^ + F^K^}\m' Jc ) , (50) 

Here 

= Y, tta(Q 2 )iWc^(p' c )Y . (51) 

In Eqs. (50), and (51) (p c - p' c ) 2 = - Q 2 , p\ = p' c 2 = 
M 2 , M c is the mass of the composite particle. In the 
state vector variables spins and masses are omitted. 

In the frame of RHD the form factors of composite 
systems ff n are to be expressed in terms of RHD wave 
functions and constituents form factors. 

In RHD a state of two particle interacting system is 
described by a vector in the direct product of two one- 
particle Hilbcrt spaces (see, e.g., Ref. [1]). So, the matrix 
element in RHD can be decomposed in the basis (30): 

/dP dP ' 
Ncg N' d ^ 

x(Pc m Jc \P ,\fs,J ,l,S, mj) 

x(P,Vi, J, Z,S,mj|i M (0)\P',V7, J', I', S', mj,) 
x (P ' , \fs', J' ,1' ,S' , mj, \p c ' , m' Jc ) . (52) 

Here the sum is over variables J ,</' ,l,V ,S ,5" ,mj ,mj> , and 
(P ' , Vs 7 ,J',l',S', mji \p' , m' Jc ) is the wave function in 
the sense of IF RHD. 
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Let us write the IF RHD wave function in the form 
slightly more general than in Ref. [1] : 

(P ,<Js,J ,l,S ,mj\p c ,m Jc ) 

= N C 6(P -Pc)5j c j5 mjcmj <pf§(k) . (53) 




The RHD wave function of constituents relative motion 
with fixed total angular momentum is defined as 

<p{§{k{s)) = ^(l-V 2 /s 2 )u ls (k) k , (54) 
and is normalized by the condition 

I uf s (k)k 2 dk = l. (55) 
is J 

Here -q = M\ — M| ,uis(k) is a model wave function. 

To calculate the form factors in Eqs. (50), and (51) 
let us write the Wigner-Eckart decomposition on the 
Poincare group for the current matrix element in the 
r.h.s. of Eq. (52) should be written. However, now there 
are some difficulties. 

In the previous sections we were dealing with the 
parametrization of local operator matrix elements in the 
case when the transformations of the state vectors and 
of the operators were defined by one and the same rep- 
resentation of the quantum mechanical Poincare group. 
While describing the composite systems in RHD a dif- 
ferent situation can arise when the state vectors and the 
operator under consideration are transformed following 
different representations of this group. 

It is just such a situation takes place in the case of 
the matrix element in the r.h.s. of Eq. (52). The op- 
erator describes the system of two interacting particles 
and transforms following the representation with Lorentz 
boosts generators depending on the interaction [1]. The 
state vectors physically describe the system of two free 
particles and present the basis of a representation with 
interaction-independent generators. So, the Wigncr 
Eckart decomposition in the form used in Section II can 
not be applied directly to the matrix element in the in- 
tegrand in the r.h.s. of Eq. (52). This is caused by the 
fact that it is impossible to construct 4-vectors describ- 
ing the matrix element transformation properties under 
the action of Lorentz boosts from the variables entering 
the state vectors (contrary to the case of, e.g., Eq. (50)). 
In fact, the possibility of matrix element representation 
in the form (50) is based on the following fact. Let us 
act by Lorentz transformation on the operator 



L>- 1 (A).^(0)?7(A)=.^(0) . (56) 
We obtain the following chain of equalities: 
{Pc,mj c \f{0)\p c ',m'j c ) 

= {p c ,mj c \U~ 1 {K)f{Q)U{K)\p c ',m' Jc ) 
= Y, (mjcliD^iR^lmjc) 

rhj c ^rh , Jc 

x (Ap c , m Jc |j"(0)|Ap c ' , m' Jc )( m' Jc \D J < (R A )\ m' Jc ). 

(57) 

Here D Jc (Ra) is rotation matrix realizing the angular 
momentum transformation under the action of Lorentz 
transformations. The equalities (57) show that the trans- 
formation properties of the current as a 4-vector (22) can 
be described using the 4-vectors of the initial and the 
final states. This means that the canonical parameter- 
ization is the realization of the Wigner-Eckart theorem 
on the Poincare group. 

In the case of the current matrix element in the r.h.s. 
of Eq. (52) the relations (57) are not valid and direct 
application of the Wigner-Eckart theorem is impossible. 

However, it can be shown that for the matrix element 
in Eq. (52) considered as a generalized function (that is 
considered as an object having sense only under integrals 
and sums in Eq. (52)), the equality (57) is valid in the 
weak sense. 

Let us consider the matrix element in question as a 
regular Lorentz covariant generalized function (see, e.g., 
Ref. [23]). Using Eq. (53) let us rewrite Eq. (52) in the 
following form: 

(Pcmjclj^lpc' ,m' Jc ) = 

/AT N' 
c c , dVSd^rf s %s)<p{, c sl (s') 
... .... N CG N CG 

x{PcVs,J c ,l,S, m Jc liju. (0) \p c ' , VP , J c , I' , S' , m' Jc ) . 

(58) 

Here it is taken into account that the current operator 
.7^(0) is diagonal in total angular momentum of the com- 
posite system. 

Let us make use of the fact that the set of the states 
(30) is complete: 
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x\P,y/s,J,l,S,mj){P,V8,J,l,S,mj\ . (59) 

Here the sum is over all the discrete variables of the basis 
(30). 

Under the integral the matrix element of the trans- 
formed current satisfies the following chain of equalities 
( (53) and (59) are taken into account): 

/N N' 
NCG N CG 

x(p c ,^,J c ,l,S, m Jc |L>- 1 (A)^(0)L>(A) 
x\p c ' ,Vs 7 , J c ,l' ,S',m' Jc ) 
= ( Vc , m Jc | U~ 1 ( A) 1 j„ (0) I U ( A) | p c ' , m' Jc ) 
= £ {mJc\[D J *(RA.)]- 1 \m Jc ) 

fhj c ,rh' Jc 

x (Ap c , m Jc \if{0)i \Kp c ' , m' Jc )( m' Jc \D J ° (R A )\ m' Jc ) 

/N N 1 

x J2 (mJc\[D J °{RA)]- 1 \m Jc ) 

mj c ,m' Jc 

x (Ap c ,y/s,J c ,l,S, m Jc \j^{Q)\Kp c ' , V7, J c , I' , S' , m Jc ) 
x(m' Jc \D J °(R A )\m' Jc ) . (60) 

It is easy to see that under the integral the current ma- 
trix element satisfies the equalities analogous to Eq. (57), 
so now it is possible to use the parameterization method 
of the previous sections under the integral, that is to use 
the Wigner-Eckart theorem in the weak sense. 

The next step is a parameterization of the matrix ele- 
ment in the r.h.s. of Eq. (58). The r.h.s. can be written 
as a functional on the space of test functions of the form 
(see Eq. (54, too)): 

i> ll ' ss '(s,s')=u l s(k(s))u l , s >(k(s')). (61) 

and Eq. (58) can be rewritten as a functional in R 2 with 
variables (s,s') (see Eqs. (47)-(48), too): 

(Pcmjclj^lpc' ,m' Jc ) 



u^s' J n cgN cg 

x{Pc,Vs,Jc,l,S, mj c \jn(0)\p c ' , V7 ,J c ,l' ,S' , m' Jc ) . 

(62) 

The measure in the integral (62) is chosen with the ac- 
count of the relativistic density of states, subject to the 
normalization (54), (55) (see Eq. (48, too)): 

d/i(s, s') = 16 0(s - (Mi + M 2 ) 2 ) 6(s' - (Mi + M 2 ) 2 ) 

x yj v^(l - v 2 /s 2 ) V^(l - rj 2 /a' 2 ) dfj,(s) dn(s') . (63) 

dfj,(s) is given by Eq. (48)). 

The sums over discrete invariant variables can be trans- 
formed into integrals by introducing the adequate delta- 
functions. Then the obtained expressions are functionals 
in R 6 . 

The functional in the r.h.s. of Eq. (62) defines a 
Lorentz covariant generalized function, generated by the 
current operator matrix element. The integral in Eq. (62) 
converges uniformly due to the definition (61). 

Taking into account Eq. (60) we decompose the matrix 
element in the r.h.s. of Eq. (62) into the set of linearly 
independent scalars entering the r.h.s. of Eq.(50): 

NcK 
N CG N' CG 

x (Pc , \fs, Jc, I, S, m Jc \j^(0)\p c ' , VP, J c , I', S', m' Jc ) 

= (m Jc \D J <(p c , p' c ) J2(iPc^(p' c )r 

xA l ^ ss '(s,Q 2 , S ')\m' Jc ), (64) 

Here A l r l lt f s (s,Q 2 ,s') is a Lorentz covariant generalized 
function. 

Taking into account Eq. (64) and comparing the r.h.s. 
of Eq. (50) with Eq. (62) we obtain: 

£ / dv( s , s i)ip u ' ss '( s , s >) 

l,l',S,S' J 

x(m Jc \A l ^ SS '(s,Q 2 ,s')\m' Jc ) 
= <m Jc |[/f n ^ + /| n rVc) 
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+fLR» + ttnK ll ]\m' Jc ) . 



(65) 



All the form factors in the r.h.s. of Eq. (65) are nonzero 
if the generalized function A contains parts that are di- 
agonal (Ai) and non-diagonal (^2) in m.jc , m [j c - 



A II' SS' 1 



(s, Q, s') = AZ%> (s, Q , s>) + A l Lf (s, Q , s>) . 



'"' S ' S "< ' > 2 s ') 
(66) 

For the diagonal part we have from Eq. (65): 
£ f d»(s,s')i; u ' ss '(s,s') 

l.l'S.S'J 



x(m Jc \A l l n l s (s,Q 2 ,s')\mj c ) 
(mjc\ [IIM K + ftM K»} \m Jc 



(67) 



The notation ff n [ip} in the r.h.s. emphasizes the fact that 
form factors of composite systems are functionals on the 
wave functions of the intrinsic motion and so, due to 
Eq. (61), on the test functions. 

Let the equality (67) be valid for any test function 



u'ss' 



(s, s'). When the test functions (the intrinsic mo- 



tion wave functions) are changed the vectors in the r.h.s. 
are not changed because according to the essence of the 
parametrization (50) they do not depend on the model 
for the particle intrinsic structure. So, when the test 
functions are varied the vector of the r.h.s. of Eq. (67) 
remains in the hyperplace defined by the vectors , K'^. 

When test functions are varied arbitrarily the vector in 
the left-hand side (l.h.s.) of Eq. (67) can take, in general, 
an arbitrary direction. So, the requirement of the validity 
of Eq. (67) in the whole space of our test functions is that 
the l.h.s. generalized function should have the form: 



Atf\s,Q\s>) = K^Gt SS \s,Q\s>) 
+K,Gt SS \s,Q\s') . 



(68) 



Here G\ l n ss (s,Q 2 ,s') , i = 1,4 are Lorentz invariant 
generalized functions. Substituting Eq. (68) in Eq. (67) 
and taking into account Eqs. (54) and (63) we obtain the 
following integral representations: 

IV S S' M i+ M 2 



x^( s )C( s ,g 2 , s ')4,( s '). 



(69) 



for i = 1,4. In the case of matrix element in Eq. (65) non- 
diagonal in mj c , m' Jc we can proceed in an analogous 



way and obtain an analogous integral representations for 

fin(Q 2 ) , * = 2,3. 

So, the matrix element in the r.h.s. of Eq. (62) con- 
sidered as Lorentz covariant generalized function can be 
written as the following decomposition of the type of 
Wigner-Eckart decomposition: 

{p c ,y/s ,J c ,l ,S , m Jc |j M (0)|p c ' , Vs 7 , J c ,1' ,S' , m' Jc ) 

= NC n!nT {mjclD ' Ic{Pc > P ' c) K + WVc) 

+TzR li +J : tK ll \\m'j c ) . (70) 



2J C 

Ti = £ Gl ss '(s, Q\ s')(ip c ^(p' c )r , (71) 

n=0 

with the constraint (69). 

In Eqs. (70), and (71) the form factors contain all the 
information about the physics of the transition described 
by the operator j^(0). They are connected with the 
composite particle form factors (50), and (51) through 
Eq. (69). In particular, physical approximations are 
formulated in our approach in terms of form factors 
G l ^ ss '{s,Q 2 ,s') (see Ref. [1] for details). The matrix 
element transformation properties are given by the 4 
vectors in the r.h.s. of Eq. (70). 

It is worth to emphasize that it is necessary to con- 
sider the composite system form factors as the functionals 
generated by the Lorentz invariant generalized functions 
Gl ss '(s,Q\s>). 

Now let us impose the additional conditions on the 
matrix elements in Eqs. (50), and (70) in the same way 
as we did in Sec. II in Eqs. (18)-(21) and in Sec. Ill in 
Eqs. (42), (44), and (45). The r.h.s. of equalities (50) 
and (70) contain the same 4-vectors and the same sets 
of Lorentz scalars (51) and (71), so, to take into account 
the additional conditions it is necessary to redefine these 
4-vectors according to Eqs. (18)— (21). For example, the 
conservation law gives T^=i). It is easy to see that for the 
redefined form factors the equality (69) remains valid. 

Let us write the parameterization (70), (71) for the 
particular case of composite particle electromagnetic cur- 
rent with quantum numbers J = J' = S = S' = 1, which 
is realized, for example, in the case of deutcron. Separat- 
ing the quadrupole form factor in analogy with Eq. (46) 
and using Eqs. (70), and (71) we obtain the following 
form: 

(p c ,Vs,Jc,l,S , m Jc \j^(0)\p c ' , Vs 7 , J c , I' , S' , m' Jc ) 



NcgK 



N C N C 



P(m Jc | J D 1 (p c ,p' c ) 



M, 



■F3 Rp, 



\m' Jc ) • 
(72) 
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fx = G%(s,Q 2 ,s') + Gf 2 ( s ,Q 2 ,s') {[ipcvVfo)] 2 



-W*- rV < )]2 }spfeJw' 



F 3 = G%(s,Q 2 ,s'). 



(73) 



We have taken into account that the equation 
G l 2i(s, Q 2 , s') — is valid in weak sense. The calcula- 
tion of the form factors in Eqs. (72), (73) in MIA will be 
discussed in the following Section. 



Our form factors in Eq. (75) can be written in terms 
of standard Sachs form factors for the system with the 
total angular momentum one. To do this let us write the 
parameterization of the electromagnetic current matrix 
element in the Brcit frame (sec, e.g., Rcf. [24]): 

(Pc,mjc\j^O)\Pc' ,m' Jc ) = G»(Q 2 ) , 
G (Q 2 )=2p c0 {(i'C)Gc(Q 2 ) 



+ 



Gq{Q 2 ) 



2M? 



(76) 



V. ELECTROMAGNETIC CURRENT MATRIX 
ELEMENT FOR COMPOSITE SYSTEMS WITH 
J =1 

In this section we make use of the formalism developed 
in the previous sections to describe composite systems 
with J = J' = S = S' =1. The case of zero total angular 
momentum and zero total spin was considered in detail 
in Ref. [1]. 

If we use the form (50) and impose the additional con- 
ditions in analogy with Eqs. (72), and (73), then the 
parameterization (50) takes the following form: 



(Pcmjclj^O^Pc' ,m' Jc ) 
= (m Jc \D\p c ,p' c ) ]T F 4 c 4|m' Jc ). 

i — 1 3 



(74) 



Here 



f\ = fUQ 2 ) +~f c AQ 2 ) {[wc^(p' c )} 2 



|sp[^rVc)] 2 } Spb ,j w) 



MQ 2 ) ■ 



(75) 



G(Q 2 



Pea 



?(CQ)-C(£'Q) G M (Q 2 ). 



Here Gc , Gq , Gm are the charge, quadrupole and 
magnetic form factors, respectively. 

The polarization vector in the Breit frame has the fol- 
lowing form: 



^(±l) = -±=(0, Tl, - i, 0) 
^(0) = (0,0,0,1). 



(77) 



The variables in £ are total angular momentum projec- 
tions. 

In the Breit frame: 



(p c - p: 



(0, Q) 



1 ^ 



(78) 



p c o = ^/M2 + -QZ, Q = (0,0,0). 

Comparing Eq. (74) with Eq. (76) and taking into ac- 
count the fact that in the Breit system D x m ; m „ (p c , p' c ) — 
S m , m " • we have: 



Al = (Pc+p'Jt, , Al = -^-e^x.p" cP ' c X r°(p' c ) . 

The form factors are redefined here (compare to Eq. (51)) 
so as to have the meaning of the charge, quadrupole and 
magnetic form factors (see Eq. (46), too). 

It should be mentioned that formally Eqs. (74), and 
(75) can be obtained from Eqs. (39) - (46) for the current 
matrix element of the free system if = p cfl7 P^ = j/ . 

The representation (74), and (75) of the matrix ele- 
ment satisfies all the conditions on the composite system 
electromagnetic current [1] (see Sec. II). 



o m 2 ~ 



G C (Q 2 ) = AW) . G Q (Q 2 ) 

G M (Q 2 ) = -M c f%(Q 2 ) 



Q 2 



(79) 



The integral representations for the composite system 
form factors in Eq. (74) are obtained in analogy with 
Eq. (69) from Eqs. (58), (72)-(75): 



/»(Q 2 ) = E / <V« dsfs 1 ip l (s) G l £ {s,Q 2 , s')</ (s' 



)• 

(80) 
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In Eq. (80) the variables S = S' =1 are omitted 

Let us pay a special attention to the following points. 

1. The chain of equalities (74)-(80) guarantees that 
the general conditions (22)-(26) for the electromagnetic 
current operator are fulfilled on each step of the calcula- 
tion. 

2. The result for the form factors (80) does not depend 
on the actual values of mj , m'j entering the l.h.s. of Eq. 
(74). In the standard light front dynamics approach such 
an ambiguity does exist [14,15]. 

3. As compared with light front dynamics [8,25] we 
have not used any particular ("good") current compo- 
nents in the calculation of form factors. 

Using the equality (79) we obtain three scalar equal- 
ities for the standard Sachs form factors of composite 
system with total angular momentum one and total spin 
one: 

G C (Q 2 ) = ]T f d^dV^ip l (s)GUs,Q 2 ,s')/(s'). 
i,v J 

Gq(Q 2 ) 

1 M 2 f 

= -n^E/ ^V^(s)Gi' 2 ( s ,Q 2 , s V (*')• 

(81) 

G M {Q 2 ) 

= -M c J2 I d^d^<p l (s)G%(s,Q 2 ,s')<p l '(s') . 
i,v 

The next case of our consideration is the case of p 
meson. In our composite quark model the u- and d- 
quarks are in the S state of the relative motion, that is 
I = I' = 0. So there is no summation in Eq. (81). 

Let us use for (81) the modified impulse approximation 
formulated in terms of form factors Gf q (s, Q 2 , s'). The 
physical meaning of this approximation is considered in 
detail in Ref. [1]. In the frame of MIA the invariant 
form factors Gf q (s, Q 2 , s') in (81) are changed by the free 
two-particle invariant form factors given by Eqs. (39), 
(46), and (A1)-(A3) and describing the electromagnetic 
properties of the system of two free particles. So, the 
equations to be used for the calculation of the p-meson 
properties in MIA are the following: 

G C (Q 2 )- / d^dV7 V { S )g nc ( S ,Q 2 ,s') V { S ') , 



2 m r 

Gq(Q 2 ) = ~qT~ J d^dV?<p(s)g Q Q(s,Q 2 ,s')<p(s') , 

(82) 

G M (Q 2 ) = -M c J d^dV?^(s)g M(s,Q 2 ,s')<p(s') . 

It should be noted that Eq. (82) can be obtained for- 
mally by changing in the r.h.s of Eq. (52) the current 
matrix element with interaction into the current of the 
free two-particle system (39), (40), (46) along with the 
following change in the covariant part (39): 

D J (P,P') , Aj, D J (P,P>) , Al | P=Pc , P , =P , , 

i = 1,2,3 (83) 

So, in our approach the current for the interacting system 
(to be precise, its matrix element) is obtained from the 
free two-particle current by including the interaction in 
the covariant part of D J , in Eq. (39). In Eq. (83) the 
components of vectors p , p' differ from vectors P , P' of 
the free two-particle system in Eqs. (39), and (40) as to 
take the interaction into account. So, our approach does 
not contradict to the general statement [8] saying that 
the current of an interacting system should be interaction 
dependent for the conditions of Lorentz covariance and 
the current conservation law be fulfilled. 

The detailed discussion of Lorentz-covariance condi- 
tion, conservation laws, cluster separability and of the 
charge nonrenormalizability, as well as of the dynamical 
content of MIA, can be found in Ref. [1]. 

VI. p - MESON ELECTROMAGNETIC 
STRUCTURE 

Electroweak properties of composite hadron systems 
were described in the RHD approach in a number of pa- 
pers. The most popular approach in the frame of RHD is 
the light-front dynamics [8,14-16,25-27]. Recently some 
calculations in the frame of instant form [28,29] and point 
form dynamics [30] have appeared. The p - meson elec- 
tromagnetic structure was calculated in Rcfs. [14,15] in 
the light-front dynamics approach. 

In this section we make use of the results obtained in 
the previous sections to calculation the p - meson elec- 
tromagnetic properties. 

The p - meson electromagnetic form factors are cal- 
culated using Eqs. (82) in MIA. The wave functions are 
taken according to Eqs. (54), and (55). We suppose that 
quarks are in the S state of relative motion. 
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For the description of the relative motion of quarks (as 
in Ref. [1] in the case of pion) in Eq. (54) the following 
phenomcnological wave functions are used: 

1. A Gaussian or harmonic oscillator (HO) wave func- 
tion (sec, e.g., Ref. [25]) 

it(fc) = N HO cxp {-k 2 /2 b 2 ) , (84) 

2. A power-law (PL) wave function [27]: 

u(jfe) = N PL (k 2 /b 2 + l)~ n , n = 2 , 3 . (85) 

3. The wave function with linear confinement from 
Ref. [31]: 

u(r) = N T cxp(-w 3/2 - fir) , (86) 



a = - sj2 M r a , (3 = M r b . 

Here a , b are the parameters of linear and Coulomb parts 
of potential respectively, M r is the reduced mass of the 
two-particle system. 

For Sachs form factors of quarks (see (A1)-(A3)) we 
have [32]: 

G q E (Q 2 ) = e q f(Q 2 ) , G q M (Q 2 ) = (e q + Kq ) f q (Q 2 ) , 

(87) 

where e q is the quark charge and K q is the quark anoma- 
lous magnetic moment. However, for f q (Q 2 ) the form of 
Ref. [33,34] is used instead of Ref. [32]: 



f q (Q 2 



l 



1 + ln(l 



(r 2 q )Q 2 /Q) 



(88) 



Here (r 2 ) is the mean square radius (MSR) of constituent 
quark. 

The details describing the cause of the choice (88) for 



the function f q (Q 2 ) can be found in Ref. [33] (see also 
Ref. [1]) and are based on the fact that this form gives 
the asymptotics of the pion form factor as Q 2 — > oo that 
coincides with the QCD asymptotics [35]. 

So, for the calculations we use a standard set of param- 
eters of CQM. The structure of the constituent quark is 
described by the following parameters: M u = Mj = M 
is the constituent quark mass, k u , Kg is the constituent 
quarks anomalous magnetic moments, ( r 2 ) = ( rV) = 
(r 2 ) is the quark MSR. The interaction of quarks in p 
meson is characterized by wave functions (84) - (86) with 
the parameters b , a , (3. 

The parameters were fixed in our calculation as fol- 
lows. We use M=0.25 GeV [36] for the quark mass. The 
quark anomalous magnetic moments enter the equations 



through the sum (k u + Kg) and we take k u + Kg = 0.09 
in natural units following [37]. 

For the quark MSR we use the relation (r 2 ) ~ 0.3/M 2 
[32,38]. 

The parameter b in the Coulomb part of the potential 
in Eq. (86) is b = 0.7867. This gives the value a s = 0.59 
for systems of light quarks [39] . 

We choose the parameters b in Eqs. (84) and (85) 
and a in Eq. (86) so as to fit the MSR of p meson. 
The p - meson MSR is given by the equation (r 2 ) — 
(rl) = 0.11±0.06 fm 2 [32,38]. For the pion MSR the 
experimental data [40] is taken: (r 2 ) 1 / 2 = 0.657±0.012 
fm. We use the following relation: 



<r 2 ) = -6G' c (0) 



(89) 



The p - meson mass in Eq. (82) is taken from Ref. [41]: 



M p = 763. Oil. 3 MeV . 



The magnetic \i p and the quadrupole Q p moments of 
p meson were calculated using the relations [24] : 



Gm(0) = ^ p , G Q (0) = M 2 p Q p 



(90) 



The static limit in Eqs. (82) gives the following relativis- 
tic expressions for moments: 



L(s) 



i r°° 

5 / dy/i 
L J2M 



+ (K U + K d ) 



AM I 



<p 2 (s) 



Vs-AM 2 



{l-L(s) 



L(s) 



(91) 



2M ^s{s~4M 2 ) 

+ K U + Kg] L{S) , 



M 



Vi + 2M 



(92) 



2M 2 



Vs-4M 2 ( v / i+2M) 



+ ln 



y/s- Vs-4M 2 



Let us note that the nonzero p - meson quadrupole mo- 
ment appears due to the relativistic effect of Wigner spin 
rotation of quarks only (see (A2)). So, measuring of the 
quadrupole moment can be a test of the relativistic in- 
variance in the confinement region. 

The Wigner rotation contributes to the magnetic mo- 
ment, too. Without spin rotation [w\ = cj 2 = in (A3)) 
we have for the magnetic moment: 
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x / dJ~s ■ , 

J 2 m y/s-AM 2 



1 - -L(s) 
2 1 ' 



(93) 



To estimate the contribution of the relativistic effects 
to the p - meson electromagnetic structure the non- 
relativistic calculation of electromagnetic form factors 
and static moments was performed. The non-relativistic 
limit of equations (82) gives the following forms of the 
form factors: 



G C (Q 2 



G E (Q 2 ) + G d E (Q 2 )) HQ 2 ) , G Q (Q 2 ) = 0, 



Gm(Q 2 ) = |f (G U M (Q 2 ) + GfAQ 2 )) I{Q 2 ) , (94) 



i(Q 2 



k 2 dk k' 2 dk' u(k) g NR (k, Q 2 , k')u(k') 



9NR(k,Q 2 ,k') = 



kk'Q 



k' 



9 k' 



Let us note that the obtained result (94) coincides with 
the one derived in standard non-relativistic calculations 
for composite system form factors in the impulse approx- 
imation [42]. So, Eq. (82) can be considered as a rela- 
tivistic generalization of the equations of Ref . [42] . 

TABLE I. The p - meson statical moments obtained with 
the different model wave functions (84) - (86). (r% R ) is non- 
relativistic MSR, (f 2 ) is relativistic MSR without spin rota- 
tion contribution, p p is relativistic magnetic moment (91), p p 
is relativistic magnetic moment without spin rotation (93), 
Q p is quadmpole moment (92). The wave functions param- 
eters are obtained from the fitting of p - meson MSR ob- 
tained through the relativistic calculation with spin rotation, 
(r 2 P ) - (rl) = 0.11±0.06 fm 2 [32,38]. The pion MSR was taken 
from the experimental data [40]: (r 2 ) 1/2 = 0.657±0.012 fm. 
The magnetic moments are in natural units, the quadmpole 
moments and MSR are in fm 2 . The parameters b in (84) and 
(85) - are in GeV and a in (86) - in GeV 2 . The values of 
other parameters are given in the text. The number of the 
significant digits is chosen so as to demonstrate the extent of 
model dependence of calculations. 



Wave 














functions 


b , a 


(r 2 NR) 


(f 2 ) 


p P 


p P 


Q P 


(84) 


0.231 


0.275 


0.731 


0.852 


0.966 


-0.0065 


(85) n=2 


0.302 


0.319 


0.711 


0.864 


0.972 


-0.0059 


(85) n=3 


0.430 


0.305 


0.710 


0.866 


0.973 


-0.0061 


(86) 


0.028 


0.301 


0.711 


0.865 


0.973 


-0.0061 



The calculation of the non-relativistic magnetic mo- 
ments using Eq. (90) gives the following result which does 
not depend on the model wave functions: 



PpNR 



1 



(95) 



The values of parameters being fixed before we obtain 
PpNR = 1-09. 

The results of the calculation of the p - meson stat- 
ical moments are given in the Table I. The relativistic 
results for the same parameters but without Wigner spin 
rotation (u>i — lo-i — in (A1)-(A3)), as well as of non- 
relativistic calculation are given, too. The number of sig- 
nificant digits is chosen so as to demonstrate the extent 
of the model dependence of the calculations. 

As one can see from the table the contributions of the 
spin rotation to the magnetic moment and to MSR de- 
pend weakly on the model for the quarks interaction in p 
meson. This contribution to the MSR is 24%-26% and is 
negative. This result differs from that of the paper [14] 
where the contribution of spin rotation (Melosh rotation) 
to MSR calculated in the frame of light-front dynamics 
is positive. 



1.00 



O 



0.10 



0.01 




Q 2 (GeV 2 ) 



FIG. 1. The results of the calculations of the p-meson 
charge form factor with different model wave functions. The 
values of the parameters are given in the Table I. The solid line 
represents the relativistic calculation with the wave function 
(84), the dashed line - with (85) for n = 3, dash-dot-line - 
with (86), dotted line - with (85) for n — 2, dot-dot-dash-linc 
- the non-relativistic calculation with (84) 
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Q 2 (GeV 2 ) 

FIG. 2. The results of the calculations of the p-meson 
quadrupole form factor with different model wave functions, 
legend as in Fig.l 
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FIG. 3. The results of the calculations of the p-meson mag- 
netic form factor with different model wave functions, legend 
as in Fig.l 
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Q 2 (GeV 2 ) 

FIG. 4. The contribution of the relativistic spin rotation ef- 
fect. The results of the relativistic calculation of the p-meson 
charge form factor with the wave function (84) using the same 
parameters as in Fig.l. Solid line represents the relativistic 
calculation with spin rotation, dotted line - the relativistic 
calculation without spin rotation. 

The spin rotation contribution to the magnetic mo- 
ment in our calculations is 11%-12% and is negative, 
too. The total relativistic corrections to MSR in our ap- 
proach are positive and enlarge the non-relativistic value 
essentially - almost twice in the case of the model (84) 
and for 70% - 80% for the models (85), (86). 

The total relativistic corrections for the magnetic mo- 
ment as compared to the non-relativistic result (see 
Eq. (95)) are negative and have the value of 21%-22%. 
Let us note, that in the light-front dynamics approach 
[15] a different result was obtained: the positive rela- 
tivistic correction of the value of 10% to the magnetic 
moment. 

Let us note especially that in our calculation of the p 
- meson quadrupole moment no ambiguities connected 
with rotation symmetry arises - in contrast with the 
light-front dynamic calculations [14,15]. 

The results of calculations for the p - meson electro- 
magnetic form factors are represented in Fig. 1-3. 

Let us note that our charge form factor has no dip in 
contrast with the results of the papers [14,15]. This prin- 
cipal difference is probably due to the effects of rotation 
symmetry breaking in Refs. [14,15]. The relativistic cor- 
rections in our approach diminish essentially the rate of 
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the decreasing of the charge and magnetic form factors 
at large values of momentum transfer. We demonstrate 
in the figures the case of the model (84) with the ex- 
ponential decreasing of the nonrclativistic form factors 
with the increasing Q 2 . The nonrelativistic quadrupole 
form factor is zero in the absence of the D-state in the 
two-particle system. 

In Fig. 4 the contribution of Wigner rotation of quark 
spins to the p - meson charge form factor is shown. This 
contribution depends weakly on the momentum transfer 
in the range from 1 to 5 GeV 2 and its value is approx- 
imately 10%. In our calculation the sign of this contri- 
bution differs from that obtained in the light-front dy- 
namics approach [14] . Let us note that similar difference 
takes place in the case of pion electromagnetic structure, 
too (see calculations in Ref. [25] and Ref. [34]). 



VII. CONCLUSION 

The method of construction of the electromagnetic 
current matrix elements for the relativistic two-particle 
composite systems with nonzero total angular momen- 
tum is developed in the frame of the instant form of RHD. 

The method makes use of the Wigner-Eckart theorem 
on the Poincare group. It enables one to extract from the 
matrix elements the reduced matrix elements - invariant 
form factors - which in the case of composite systems are 
generalized functions. 

The obtained current operator matrix elements satisfy 
the Lorentz-covariance condition and the conservation 
law. 

The modified impulse approximation — with the phys- 
ical content of the relativistic impulse approximation — 
is formulated in terms of reduced matrix elements. MIA 
conserves Lorentz covariance of electromagnetic current 
and the current conservation law. 

The developed formalism is used to obtain a reasonable 
description of the static moments and the electromag- 
netic form factors of p meson. A number of relativistic 
effects are obtained, for example, the nonzero quadrupole 
moment (in the case of S state) due to the relativistic 
Wigner spin rotation. 

So, in conclusion, it is shown that the instant form 
of RHD can be used to obtain an adequate description 
of the electroweak properties of composite systems with 
nonzero total angular momentum. 
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Appendix 1 

The charge form factor for free two-particle system: 



9oc(s,Q 2 ,s')^-R(s,Q 2 , S ')Q 2 



x{(s + s' + Q 2 )[g u e (Q 2 ) + G%(Q 2 ) 
x [2 cos(o>! — oj 2 ) + cos(cj! + UJ 2 )] 



-— £(*,qV) [giaq^ + gUQ 2 ) 

x [2 sin(cj! — u 2 ) — sin(cj! + uo 2 )]} , (^41) 
The quadrupole form factor for free two-particle system: 



1 



9 oq(s,Q 2 , S ') = ^R(s,Q 2 , S ')Q 2 



x[(s + s' + Q 2 ) [g u e (Q 2 ) + G%{Q 2 ) 



x [cos(wi — u 2 ) — cos(a>i + L0 2 )\ 



>V) 



GlAQ^ + GUQ 2 ) 



x [sin(wi - u 2 ) + sin(wi + lu 2 )]} , (A2) 
The magnetic form factor for free two-particle system: 
9om{s,Q 2 ,s') = -2R(s,Q 2 ,s') 

x fa, Q 2 , a') [gKQ 2 ) + G%{Q 2 )] sin( Wl - w 2 ) 



+ 



AM 



G U M {Q 2 ) + G d M {Q 2 ) [{s + s' + Q 2 )Q 



'3 1 
x ( - cos(wi - uj 2 ) + - cos(a>i + u 2 ) 



— A ^,Q 2 ,s') 



(V7+ 2 M)(s -s' + Q 2 ) + {s'-s + Q 2 )VJ 



+ 



(y/I + 2 M)(s' - s + Q 2 ) + (s - s' + Q 2 )yfs 



Vi(Vi + 2M) 



x [sin(o>i — cj 2 ) — sin(wi + u 2 )] 



:£ 2 (*,QV) 



v / s 7 (v / s 7 +2M) yfs{yfs + 2M) 
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x [cos(wi - uu 2 ) + cos(a>i + 0J2)]]} ■ (A3) 



Here 



R(s,Q 2 ,s') = 



(s + s' + Q 2 ) 



2y/(s -4M 2 )(s'-4M 2 ) 



0(*,QV) 



[A( S ,-Q 2 , S ')] 3/2 v 71 + Q 2 /4M2 



£(s,QV) = v/ SS 'Q 2 -M 2 A(s,-Q 2 ,s') 
Wi t^2 ~ the Wigner rotation parameters: 

£(*,QV) 



Wi = arctan ■ 



cl>2 = arctan 



a( S , S ')£( S ,Q 2 ,^ 



M(s + s> + Q 2 )a(s, s') + V^(4M 2 + Q 2 ) ' 



here a(s, s') = 2M + ^+ \/7, i?(s, Q 2 , s') = 0(s' - si) 
0(s' — S2), - the step-function. 



si,2 = 2M 2 



2M 2 



(2M 2 + Q 2 )(s - 2M 2 ) 



T^VQHQ 2 + 4M 2 ) S ( S - 4M 2 ) . 

M - the mass of u— and d quarks. The functions 
s i,2(s, Q 2 ) give the kincmatically available region in the 

plane (s, s'). M (Q 2 )- Sachs form factors of u- and d 
quarks. 
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